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Abstract 

Recovering a function from its spherical Radon transform with centers of spheres of integration restricted 
to a hypersurface is at the heart of several modern imaging technologies, including SAR, ultrasound imag¬ 
ing, and photo- and thermoacoustic tomography. In this paper we study an inversion of the spherical Radon 
transform with centers of integration restricted to cylindrical surfaces of the form TxR m , where T is a hyper¬ 
surface in R". We show that this transform can be decomposed into two lower dimensional spherical Radon 
transforms, one with centers on T and one with a planar center-set in lft" ! 11 . Together with explicit inversion 
formulas for the spherical Radon transform with a planar center-set and existing algorithms for inverting the 
spherical Radon transform with a center-set T, this yields reconstruction procedures for general cylindrical 
domains. In the special case of spherical or elliptical cylinders we obtain novel explicit inversion formulas. 
For three spatial dimensions, these inversion formulas can be implemented efficiently by backprojection 
type algorithms only requiring (9(N 4 ^ 3 ) floating point operations, where N is the total number of unknowns 
to be recovered. We present numerical results demonstrating the efficiency of the derived algorithms. 
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1. Introduction 

Let T be a hypersurface in Ift". In this paper we study the spherical Radon transform with a center-set 
T x R m that maps a function / : R' ,+m —» R to the spherical integrals 

(M f)( x , y , r ) :=-L- f f((x,y) + rco) dS(rn) for (x,y, r) e T x R" i+1 . 

|S" +m I Jgn+m-l 

Here \r\ and (x,y) e T x R" 1 are the radius and the center of the sphere of integration, respectively, §' !+m_1 
is the unit sphere in R' ,+m , and |§' ,+m ~ 1 | is the total surface area of §" +m_1 . Note that x,y as subscripts 
of (M X ,yf){x, y, r) indicate the variables in which the spherical Radon transform is applied, while in the 
argument they are placeholders for the actual data points. For example, evaluating the spherical Radon 
transform at ( x,y , r) — (0,0,0) gives (M A - v /)(0,0,0). Similar notions will be used for auxiliary transforms 
introduced below; these suggestive notations are used to facilitate the readability of the manuscript. 
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Recovering a function from its spherical Radon transform with centers restricted to a hypersurface is 
crucial for the recently developed thermoacoustic and photoacoustic tomography [[li 2J|. It is also relevant 
for other imaging technologies such SAR imaging [3|,|41] or ultrasound tomography [50. 

Explicit inversion formulas for reconstructing / from its spherical Radon transform are of theoretical 
as well of practical importance. For example, they serve as theoretical basis of backprojection-type recon¬ 
struction algorithms frequently used in practice. However, explicit inversion formulas are only known for 
some special center-sets. Such formulas exist for the case where the center-set is a hyperplane 0-11] or a 
sphere il2l4l6ll. More recently, closed-form inversions have also been derived for the cases of elliptically 
shaped center-sets (see [ 17U22I1). certain quadrics [23|,[24], oscillatory algebraic sets 112511 . and corner-like 
domains 1260 . 



Figure 1.1: Left: A phantom consisting of a superposition of three balls. Right: 3D reconstruction from its spherical Radon 
transform on a cylinder in R 3 using the inversion formula 10 derived in Section[3] 


1.1. Main contribution 


In this paper we present a general approach for deriving reconstruction algorithms and inversion formu¬ 
las for the spherical Radon transform on cylindrical surfaces yielding inversion formulas for the center-set 
T x W n , provided that an inversion formula is known for the center-set T. Our approach is based on the 
observation that the spherical Radon transform with a center-set T x R'" can be written as the composition 
of a spherical Radon transform with a center-set TcM" and another spherical Radon transform with a pla¬ 
nar center-set in R m+1 (see Theorem 12.31) . Recall that inversion formulas for the spherical Radon transform 
with a planar center-set are well known. Consequently, if an inversion formula is available for the center-set 
T c R", then this factorization yields an inversion formula for M XJI . As explicit inversion formulas are 
in particular known for spherical and elliptical center-sets, we obtain new analytic inversion formulas for 
spherical or elliptical cylinders (Theorems 13. II and 13.21) . A reconstruction result with one of our inversion 
formulas for an elliptical cylinder is shown in Figure 11.11 

We emphasize, that our factorization approach can also be combined with any other reconstruction 
algorithm for the spherical means, such as time reversal 127142911 or series expansion methods [30j,[31[], that 
can be used for the case where T is the boundary of a general bounded domain. 
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1.2. Outline 

The rest of this paper is organized as follows. In Section [2] we show that the spherical Radon transform 
with centers on T x R m can be decomposed in two partial spherical Radon transforms, one with a center- 
set T and one with a planar center-set in M m+I ; see Theorem 12.31 In Section [3] we apply the factorization 
approach to derive novel explicit inversion formulas for the case where T x 1ft."' is an elliptical or circular 
cylinder; see Theorems 13.11 and 13.21 In Section |4] we derive filtered backprojection algorithms based on the 
inversion formulas for elliptical cylinders in R 3 and present numerical results. The paper concludes with a 
short discussion in Section [5] 

2. Decomposition of the spherical Radon transform 

We study the spherical Radon transform with centers on a cylinder T x R m , where T c R" is any smooth 
hypersurface with n, m 6 N and n > 2. It will be convenient to identify R' !+m - R" x R'" and to write any 
point in R" +m in the form (x,y) e R” x R m . 

Definition 2.1 (Spherical Radon transform with a cylindrical center-set). Let / e L* (R" +m ). The spherical 
Radon transform of / with a cylindrical center-set T x R'" is defined by 

M XJ f: T x R m+1 -» R 

(x,y,r) h-> —-— I f((x,y) + rm)dS(m). 

In the following we show that M x>y can be written as the product of two lower dimensional spherical 
Radon transforms, one with a center-set T in R" and one with a planar center-set in R m+1 . These partial 
spherical Radon transforms are defined as follows. 

Definition 2.2 (Partial spherical Radon transforms). 

(a) For / e L 1 1 oc (R' ,+m ) we define 

M x f: TxR m+1 -> R: (x,y,s) ^ f f(x+soj 1 ,y)dS(a> 1 ). 

|s> I Js"-‘ 

(b) For li e L' ()c (r x R"' +l ) we define 

: T X R m+1 -» R: (x, y, r ) f h(x, (y, 0) + rcr) dS(cr). 

|S I Js m 

The operators M r and M v _ v are partial spherical Radon transforms, where the subscripts indicate the 
variables in which the integration is applied. By definition, M XJ f, M x f, and M. y s h are even in the last 
variable. 

Theorem 2.3 (Decomposition). For every f e L 1 1 oc (R” + ' ,! ) we have 

1 IS m l IS" -1 1 

MxJ = C tun \r\ l -''M > , s (\s\ n - l M x f) with C njn := - ■ (D 
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Proof. We use standard spherical coordinates <b: (0, In') x (0, n) n 2 x (0, n) m —> §' ,+m 1 written in the form 
<D(ar,/3) = (0| (o-) sin(/i|)... sin (fi m ), 4> 2 (/3)) for n e (0, 2 tt) x (0, and /3 e (0, 7 r) m , where 



sin(ai) sin(cr 2 ) • • • sin(o'„_i) 


cos(/3i) sin(/3 2 ) sin(/3 3 ) • • • sin(/3 m ) 


cos(ai) sin(ar 2 ) • • • sin(a„_i) 


cos (J3 2 ) sin(/3 3 ) • • • sin(/3 m ) 

®i(a) := 

cos(a„_ 2 ) sin(a„_i) 

and ® 2 (/3) := 

cos(/3 m _i) sin(/3 m ) 


cos(a„_i) 


cos (p m ) _ 


Notice that <I>| : (0, 2n) x (O.n)"' 2 —> §" _1 is the standard parameterization of §'' -1 using spherical co¬ 
ordinates. Further, not that dS ;! _i(a) = sin(a 2 ) ■ ■ • sin(a„_i ) n_2 do- is the surface element on §" 1 and 
sin(/ii )" _1 .. .sin(yS m )” +m_ 2 dS„_i(a)djS the surface element on gn+m i Expressing the spherical Radon 
transform M >vs - in terms of integrals over the parameter set (0, 2n) x (0, n) n ~ 2 x (0, n) m therefore yields 


(M X , y f)(x,y,r) 

+m_1 1 j(0,n) m J(0 , 


|gw+m- 

IS ”" 1 

|g«+m- 

|S" —1 1 |r| 1 " 


2n)x(0.n)” 


f{x + r sin(/3i)... sin(/3 m )<f>i(a),;y + rd> 2 (/3)) 


x sin(/3i)" 1 ... sin {/3 m )’ l+m 2 dS„_i(a) d/3 


|§ n+m 

|S"- 1 ||r| , - n 


L f 

J | J (0,7r) m 
1 —n n 

1 I J( 0,JT) m 


(M x f) (. x,y + r® 2 (/3), r sin(/3i) • • • sin(/3 m )) sin(/3i ) n 1 ... sin(/3 m )' !+m 2 d/3 


(M. r /) (x,y + r® 2 (/3), r sin(/3i)... sin(/3,„)) \r\ n 1 sin(/3i) M 1 ■ • • sin(/3 m )' 


\n—l 


\n+m —2 


dp 


l§" 


— f (k|" ! M r /) (x,y + r® 2 (/3), r sin(/3i)... sin(/3 m )) sin(/3 2 )... sin(/3 m ) m 1 
I J( 0,n) m 


dp. 


Next notice that p i-> (<I> 2 (/3), sin(/3i)... sin(/3 m )) is a parameterization of S+ := § m n {(y, 5 ) e M m+1 : s > 0) 
and that sin(/3 2 )... sin (J3 m ) m ~ l dp is the corresponding surface element. Because 1 M v / is even in s, 


(M x. y f)(x,y, r) 

\ r \l- n |gfi-l| 




l§"™ _± l J(0,2?r)X(0,7r) m_1 
. 1 1 —n ii&mi ii&n-li 


OI” *M X f)(x,y + r<h 2 (/3), r sin(/3i) • • • sin(/3,„)) sin/3 2 • • • sin(/3 m ) m 1 dp 


which is the desired identity. 


□ 


According to Theorem 12.31 a possible strategy to recover the function /: R" + "' —> R from its spherical 
Radon transform M XJ f is by first inverting M_ VjJ and subsequently inverting M v . Formulas and algorithms 
for inverting the spherical Radon transform M y>i with centers on a hyperplane are well known; see tun. 
If T is the boundary of a bounded domain B c M" ! , then efficient reconstruction algorithms for recovering 
a function with support in B are known, including time reversal [25, 27-2^] and series expansion methods 
sm. Hence in such situations Theorem 12. 3 1 yields efficient reconstruction algorithms. 

Further, explicit formulas for inverting M v are known for some particular surfaces T, including spheres 
ifMl and ellipsoids 1 17-22]. Combining inversion formulas for M,_ v and M T gives us explicit formulas 
for recovering a function from its spherical Radon transform on elliptical cylinders. In the following section 
we derive such formulas for elliptical cylinders. 
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3. Explicit inversion formulas for elliptical cylinders 

In this section we consider the case where T = dE,\ is the boundary of the solid ellipsoid E A '■= {x e 
R" : |A“ 1 x| < 1), where A is a diagonal n x n matrix with positive diagonal entries a\,... ,a n . We use 
the factorization of the spherical Radon transform to derive explicit inversion formulas for elliptical and 
spherical cylinders. 

We define the partial spherical back-projections of h: dE,\ x R m+1 —> R by 

(M l tS h)(x,y, s ) := J" h |x, y', yj\y — v'| 2 + \ 2 j dv' (x,y, s ) e dE A x R m+I , 

(M| h)(x,y) := f h(Acr,y, \x - Aoj)dS(cr) (x,y) e E A ~X R m - 

JS "- 1 


Note that the integral in the definition of M \ s h may not converge absolutely. In such a situation M( V( h can 
be defined via extension to the space of temperate distributions; see |@] for details. Further, we denote by 
A V i the Laplacian with respect to (y,s), A Ax Z"= i(lthe Laplacian with respect to Ax, H v the 
Hilbert transform, and by D v := (2s)~ ] d s the derivative with respect to .v 2 . 

Theorem 12.31 in combination with the inversion formulas for the spherical Radon transform on a pla¬ 
nar center-set and an elliptical center-set dE A yields explicit inversion formulas for For example, 

combining the inversion formulas of [@] and [19, 22], we get the following result. 


Theorem 3.1 (Inversion formula for elliptic cylinders). For f e C‘^(E A x R m ) and (x,y) e E,\ x R m , we 
have 

2 " 2 dCt(A ' |S " + "^ \r\ n ~ l M x , y f) (x, y ), 


f(x,y) = 


(;n - 2)\{2n) m |S»- 2 ||§» 


with 


f(-1)”/ 2 (x, v,s) ifn > 2 is even, 

(B S h)(x,y,s):=l . _\ n ( \ 

|(_ 1 )(« i)/2 ^ s p)« 2 H s h) (x,y, s ) ifn > 3 is odd. 


Proof. See Appendix A.l 


( 2 ) 

(3) 

□ 


In the case of a circular cylinder with n = 2 we have the following simpler inversion formula. 

Theorem 3.2 (Inversion formula for circular cylinder). Let Dr c R 2 denote a disc of radius R centered at 
the origin and let f 6 C“(D^ X R m ). Then, for every (x, y) e Dr X R'", we have 


f(x, y) = - 


IS' 


m +11 


2(2 tt)" 


- |r| M x j) (x, y ). 


(4) 


Proof. See Appendix A.2 


□ 


In the case that m > 1 is an odd number, then (-A vvs ) ( " ! ~ l) ’ /2 is a differential operator and therefore 
local. Consequently, in such a situation, the inversion formula © is local: Recovering the function / at 
any point (x, y) uses only integrals over those spheres which pass through an arbitrarily small neighborhood 
of that point. In the case that m > 2 is an even number, then (-A v v ) (m_l ^ 2 is a non-integer power of the 
Laplacian. In such a situation, (HJl is non-local. The same argumentation shows that for even n. the inversion 
formula (0 is also local for m odd and non-local for m even. Such a different behavior for even and odd 
dimensions is also well known for the standard inversion formulas of the classical Radon transform, as well 
as for known inversion formulas for the spherical Radon transform with centers on spheres or ellipses. 
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Remark 3.3 (Elliptical cylinders in R 3 ). Consider the important special case of an elliptical cylinder 
Ea x R in R 3 , where Ea {(xi,;t 2 ) £ R 2 : (xi/ci \) 2 + (x 2 /fl 2 ) 2 < 1}- Then (0 simplifies to / = 
- det (A) /{2n) A A _ Y M Y Mf v5 |r| M Y> , which can be written in the form 


f(x,y) = J^Q r l M x,yf)\ Ao ',y', -yj\y - y '| 2 + 


dy' 


dS(cr). 


(5) 


s=|;t— Acr\ 


In Section|4]we present numerical results using ©. There we also numerically compare © to the universal 
backprojection formula (UBP) 


f(x,y) = ^ J aE v * ’ (* “ x> ) Jf ( r ld r r 1 d r rM Xt yf) {x',y, y]\y-y '\ 2 


+ s 2 dy' 


dSOc'), (6) 


- 5=p-.v'| 


where v Y is the outer unit normal to dE A . In lES, the UBP has been shown to provide an exact reconstruction 
for spherical means on planes, spheres, and circular cylinders in IP' 3 . In 114, Ts|, 2(], 23, 24], this result has 
been generalized to spheres, elliptical cylinders, and other quadrics in arbitrary spatial dimensions. ♦♦♦ 


4. Numerical implementation 

In the previous section we derived several inversion formulas for recovering a function supported in 
a cylinder. In this section we report the numerical implementations of these formulas and present some 
numerical results. We restrict ourselves to the special important case of elliptical cylinders E A x R in R 3 . 


4.1. Discrete reconstruction algorithm 

In the numerical implementation the centers have to be restricted to a finite number of discrete samples 
on the truncated cylinder 8E A x \—H, H] of finite height. Further the radii are restricted to discrete values 
in [0, /•() ]. Here 2 H is the height of the detection cylinder and the maximal radius i A ) is chosen in such a way 
that (M Xt yf)(x,y, r) = 0 for (x,y) e E A x [-H, H] and r > tq. We model these discrete data by 


g[k, m, 1] := (M XJ f) x[k], H 



for (k, m, 1) 6 {0,..., K - 1} X j-L__ L) X {0,..., M}, 


where x[k] := {a\ cos (2/rk/K), a A sin (2/rk/K)) are the angular detector positions. The aim is to derive a 
discrete algorithm, that outputs an approximation 


f[ni,n 2 ,n 3 ] 


,, n, n 2 n 3 

N t Nv l 


for (ni,n 2 ,n 3 ) e {—N Y ,..., N v }“ x {—L/2,..., L/2|. 


We implemented the inversion formulas © and © by replacing all filtration operations and all partial 
backprojections by a discrete approximation. We shall outline this for the inversion formula ©, which 
reads / = - det (A) / (27r)( A M Y M Y _ A rM X jf ) . For the numerical implementation, the operators A^*, m|, 
M yvs , and r are replaced by finite-dimensional approximations and M XJ f is replaced by the discrete data g 
resulting f := - det (A) /(2/r)(A A M Y My S r)g. Here r is a discrete multiplication operator, A A is a finite dif¬ 
ference approximation of A Ay , and and My s are obtained by discretizing M Y and M Y , v with the composite 
trapezoidal rule and a linear interpolation as described in [ 12i.l321. 

The partial backprojection operators are two-dimensional discrete backprojections for fixed k and n 3 , 
respectively. Assuming K, L, M = <9(N V ), both of these back-projection operators can be implemented using 
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Figure 4. 1: Log-log plot of the actual reconstruction time T(N) of our Matlab implementation using i[5]i in dependence of the number 
of unknowns for N v = 100,125,..., 300. A comparison with the slopes of the log-log plots of C * N 4/3 and D * N 5/3 demonstrates 
the theoretical behavior T(N) = D(N 4/3 ). The constants C, D have been selected in simply in such a manner, that the above curves 
do not interfere with each other. 


<9(N 3 ) floating point operations (see [12, 32]). Therefore, the numerical effort of evaluating ([5]) or ([6]) is 
(9(N 4/3 ), where N = (2N V + 1) 2 (L + 1) is the total number of unknowns. Note that the direct implementation 
of a standard 3D spherical back-projection operator requires (9(N 5/3 ) floating point operations, see 1132143411. 
Typically we have N v > 100, and therefore our implementations are faster than the standard 3D spherical 
back-projections by at least two orders of magnitude. Figure |4~T1 shows a log-log plot of the reconstruction 
time using ([5]) for N v e [100,125,..., 300} (where we have chosen K, L, and M proportional to NJ verifying 
that the number of floating point operations is 0(N 4 ^ 3 ). For comparison purpose, Figure |4~T1 also shows 
log-log plots of N C * N 4/3 and N i-» D * N 5/3 for some constants C, D. 


4.2. Numerical results 

For the numerical simulations presented below we consider superpositions of radially symmetric in¬ 
dicator functions. Figure 14.21 shows the considered phantom and the numerical reconstruction using the 
inversion formulas ([5]) and © for a\ = 1, a 2 = 0.8, H = 2 and ro = 4. For the discretization we used 
K = 256, L = 200, M = 400, and N v = 100. In particular - , the step size H/L = 0.01 in the vertical direction 
is chosen equal to the step size aj/N* = 0.01 in the horizontal direction and also to the radial step size 
ro/M = 0.01. We thereby use approximately five times more data points than unknowns, which mainly 
accounts for the fact that the detector elements cover a larger area than the reconstruction volume; while 
the step sizes are the same. Reducing the number of data points would introduce artifacts due to spatial un¬ 
dersampling; see [35]. One can see that the reconstruction results for all inversion formulas are very good, 
especially in the horizontal direction. In the vertical slice the vertical boundaries of the reconstructed balls 
are blurred. Such artifacts are expected and arise from truncating the observation surface, see fl, EM. 
All computations have been performed in Matlab. The three dimensional reconstructed data sets consisting 
of (2N 2 + 1)(L + 1) = 8 120 601 unknowns and K(2L + 1)(M + 1) = 41 165 056 data points using either © or 
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Figure 4.2: Reconstruction results using simulated data. Top row (from left to right): Horizontal slices of the phantom, the 
reconstruction using and the reconstruction with the UBP J6). Bottom row: The same for the vertical slices. 


© have been computed in 6 minutes on a MacBook Pro with 2.3 GHz Intel Core i7 processor. 



Figure 4.3: Reconstruction results using noisy data. From left to right: Horizontal slice using 0. horizontal slice using <D, 
vertical slice using and vertical slice using <H). 

In order to illustrate the stability of the derived discrete back-projection algorithms with respect to the 
noise we applied all algorithms to simulated data, where we added Gaussian white noise with variance equal 
to 2% of the maximal absolute value of g. The reconstruction results (using a\ - 1, a 2 - 0.8, ro = 4, H = 2, 
K = 256, L = 200, M = 400, and N v = 100 as above) for noisy data are shown in Figure 1431 As can be seen 
the implementations of all back-projection formulas perform quite stably with respect to noise. The slight 
amplification of noise is expected due to the ill-posedness of the inversion of the spherical Radon transform 
reflected by the two derivatives in any of the inversion formulas. The sensitivity with respect to noise could 
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easily be further reduced by applying a regularization strategy similar to SS3] for the spherical Radon 
transform with centers on a sphere. 

5. Conclusion 

In this paper we studied an inversion of the spherical Radon transform in the case where the centers 
of the spheres of integration are located on a cylindrical surface T x M m . We showed that this particular 
instance of the spherical Radon transform can be decomposed into two lower dimensional partial spherical 
Radon transforms, one with a center-set T c R" and one with a planar center-set in M m+I , see Theorem 
12.31 This factorization was used to derive explicit inversion formulas for elliptical and circular cylinders. 
We emphasize, that our factorization approach can also be used for more general cylinders in combination 
with existing algorithms for an inversion from spherical means on bounded domains, such as time reversal 
or series expansion. 


Appendix A. Proof of explicit inversion formulas 

Appendix A.l. Proof of Theorem 13. 1 1 

In this appendix we derive the inversion formula ([2]) stated in Theorem 13.11 For that purp ose we first 
give an inversion formula for spherical means on an ellipsoid that essentially follows from [19, 22]. 

Lemma Appendix A.l. For f e C“(£a x M m ) and every (x,y) € E A x W\ we have 


jn -3 


f(x,y) = 


det(A) 


|§"- 2 |(« - 2)! 


(Aa.tMJc.sMv/) (x, y) , 


(A.l) 


with 


| f 0 °° s'h(x,y, s')log |s 2 - s' 2 |ds' n - 2, 


(C s h) (x,y, s ) := - 


(_!)(« 1 )/ 2 ( 5 D" 2 s n 2 li)(x,y,s) n>3odd, 

(—!)("— 2 )/2 (h s D"-V- 2 /z) (x,y, s ) n> 4 even . 

Proof. For n = 2 and n = 3. the inversion formula (IA.1I) has been derived in |22] and for odd 72 > 3 in Hill . 
For even n > 4, in 119(1 the following inversion formula has been shown: 


/<«) = 


f - f> 

J S n~l n Jo v 


i «-2 g n ~2 


M x f) ( Acr,y , s)log ||jc - Acr| 2 - s 2 | sdsdS(cr). 

(A.2) 


\S n ~ 2 \(n - 2 )! 

Using 2 s(s 2 - | Act - v| 2 )' 1 = (s — \Act - x\)~ l + (s + \Act - x|) _1 , the inner integral in (1A.2I ) can be written as 

2 f°° 

(d; -2 s" - 2 M v /) (Acr,y, s) log ||jc - Acr| 2 - s 2 | s ds 

^ Jo 


{Acr,y, s ) 


i r°° 

= - P-V. I (d”“ 3 / !_2 Mt/) 

X J 0 v 7 

>.v.f 

Jo 

= -P.V. f 

n Jr s-\x- Act\ 

= (H v DrV- 2 M. v /) (Ao~,y, \x - Atr|) . 


2s ds 


s 2 - \x - Act | 2 


(D”- 3 s” _2 Mt/) {Act , y, s) 
s - \x - Aa\ 

(D"“ 3 s"- 2 M,/ ) (Act, y , s) 


ds + 


r 

Jo 


(D n - 3 s n - 2 M x f)(Acr,y,s) 


s + \x - Acr| 


ds 


ds 


Inserting the last displayed equation in (1A.2I) shows (lA.lh for the case of even n > 4. 


□ 
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Lemma Appendix A.2. For even and differentiable h: R->Kw have H V D v H v /j = -D s h. 


Proof Since h is even, H v /; and D v H v /z are odd in s. Now, for any odd function tp : 

(I 


£, we get 


«**) = i p v. r *2 *■=i (p.v. r *' -p.v. r 

n Jr s-s' n\ Jo s - s' J 0 




S + s' 


i r° 

-p.v. 

n Jo 


2 s' 


s- - (s') 


—f(s)ds . (A.3) 


Applying (1A.3I) once with <p = D. S H s h and once with <p = H v /; yields 

1 r°° 2 s' 

(H s D,HJi)(5) = - P.V. -z -—(D. s HJi)(/)ds' 

7T Jo s 1 - (s') 1 

2 s' 


i r 

Di — P.V. 

K Jo 


S 2 - (S') 


—^(H s h)(s )ds = (DjH^H S h)(s) = -D S h(s). 


□ 


Now we are ready to prove Theorem 13. II For that purpose, let / e C‘f(E A x M'"). From Theorem 12.31 
and the inversion formula [(dj] for the spherical Radon transform with a planar center-set (see also [9;, Qoj|), 
we obtain 

|§/i+m-li 


MJ = 

Together with (IA.1I) this gives 
2”~ 3 det(A) 


(27T) m |S' 


(-A y J m - l)/2 U s dMl M"' 1 M *,y/) ■ 


(A.4) 


f(x,y) = 


|S"~ 2 |(« - 2)! 

2” -3 det(A) IS"*'"- 1 


|S”- 2 |(/7 — 2)! (27r) m |S‘ 
2' !-2 det(A) |S" +m 1 | 


(a Ax mIc s M X f)(x,y) 

^7 (Aa,-M5.C, \s\ l ~ n (-A, ; ,) (m - 1)/2 H,J S M^ \r\ n ~ l M X J) (x,y) 


(2n) m (n - 2)\ |§"- 2 ||§‘ 


^ (Aa,M#B s (-A y J m -^ 2 Ml |r|” _1 M x j) (x,y ), 


with 


(B s h)(x,y, s) = 


\ (“(d.vHJi) (x,y, s') log |s 2 - V 2 | ds' if n = 2, 
i(_l)(«-P/2 (.sD'r^v-'^HJi) (x,y, s) if n > 3 odd, 
i(_l)(«-2)/2 (H 4 .D;r 3 s~ l d s H s h) (x,y, s) if n > 4 even. 


It remains to verify, that B can be written in the form ([3]). For n > 3 odd this follows from D v = 1/2 s 1 d s . 
In the case n = 2 integration by parts and applying (IA.3I) shows 

1 r°° 2 s' 

(B s h) (x,y, s) - — P.V. (H s h) (x, y, V) -r-- ds' = -(H,H,7i) (x, y, s) = h (x, y, s) . 

n Jo s 2 - s’ 2 

Finally, for the case that n > 4, repeated application of Lemma Appendix A.2 implies that B is given by 
©, and concludes the proof. 

Appendix A.2. Proof of Theorem 13. 2 1 

To prove Theorem l3.2l we use the following modification of one of the formulas of [ 12, Corollary 1.2]. 
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Lemma Appendix A.3. For f e C“(D# x M m ) and (x,y) e Dr x M m , we have 


f(x,y ) = \ (M#HA \s\ M x f) (x,y). 


(A.5) 


Proof. The proof is based on [ 12] and a range condition for M r derived in [42]. The range condition for the 
spherical Radon transform of [42] yields 


0 = 


2s 


s 2 - \x' - x\ 2 


M X f{x',y, s) ds dS(x') 


f pv T 

JdD R JO 

f fujw.y. 

Jod r Jo s -\x - x\ J 9Dr Jo 

f P.V. f sign(s)M x f(x',y,s) 

Jsd r Jr 


s) 


ds dS(x') 
s + \x’ - x\ 


ds dS(x') 
s -\x' - x\ 


= (m|h v sign(i)M A ./) (x,y). 

By the product rule, we have 

(m}hA \s\ M,/) (x,y) = (M*H, sign(.v)M,/) (x,y) + (M'H S \ s \ d s M x f) (x,y ). 


(A.6) 


According to (IA.6I) . the first term in the last displayed equation equals zero, and according to [ 12, Eq. (1.6)] 
the second term is equal to 2f(x,y). This completes the proof. □ 

The inversion formula <(4]) now follows from (I A. 41) and (IA.5I) . In fact, 


f(x,y) = ^(Mln s d s \s\M x f)(x,y) 

= ” 20 ^ ^ M *,yf) ( a >0 , 

where we used two identities H s <9v(—A^,,)-'" -1 ^ 2 = (-A ViS ) (m_1 ^ 2 J i H i and H v H s /t = -h. 
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